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SUMMARY INTRODUCTION 


A review and statistical analysis 
of the ultrasonic velocity method for 
estimating the porosity fraction in 
polycrystalline materials is presented. 
Initially, a semi-empirical model is 
developed showing the origin of the 
linear relationship between ultrasonic 
velocity and porosity fraction. Then, 
from a compilation of data produced by 
many researchers, scatter plots of 
velocity versus percent porosity data 
are shown for AI 2 O 3 , CuO, MgO, porcelain- 
based ceramics, PZT, SiC, S 1 3 N 4 , steel, 
tungsten, UO 2 , (Uo. 30 Pu 0 . 70^, and 
YBa 2 Cu 307 _ x . Linear regression analy- 
sis produced predicted slope, inter- 
cept, correlation coefficient, level 
of significance, and confidence inter- 
val statistics for the data. Velocity 
values predicted from regression analy- 
sis for fully-dense materials are in 
good agreement with those calculated 
from elastic properties. 


The physical behavior of compo- 
nents manufactured from polycrystalline 
materials is in many cases directly 
dependent on the porosity fraction 
(volume fraction of pores). As exam- 
ples concerning key properties of 
technologi cal ly- Important materials, 
porosity fraction has been shown to 
affect: ( 1 ) the strength, toughness 

and modulus of structural and refrac- 
tory materials such as steel (ref. 1 ), 
tungsten (ref. 2), SiC (ref. 3), S i 3 N 4 
(ref. 3), and AI 2 O 3 (ref. 3), (2) the 
strength of nuclear fuel materials 
such as UO 2 (refs. 4 and 5), (3) the 
thermal shock behavior and strength of 
porcelain-based ceramics (refs. 6 and 
7), (4) the dielectric and elastic 
properties of piezoelectric materials 
such as PZT (ref. 8 ), and (5) the crit- 
ical current density, diamagnetic 
response, and modulus of superconduct- 
ing ceramics such as YBa 2 Cu 307 _ x 
(refs. 9 to 11). In the latter case, 
reference 9 has shown that porosity 
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fraction variations on the order of 
1 percent in YBa 2 Cu 307 _ x samples can 
result In an order of magnitude varia- 
tion in critical current density. In 
such cases where physical properties 
are directly dependent on porosity 
fraction, the measurement of porosity 
fraction becomes important in the qual- 
ity assurance process for the material. 

Currently, various methods are 
available for measuring the porosity 
fraction of polycrystalline materials. 
The most common include dry-weight 
dimensional and liquid immersion 
(ref. 12). Other methods for obtaining 
porosity fraction include estimates 
from optical areal analysis measure- 
ments (ref. 13) and estimates from 
x-ray attenuation measurements 
(ref. 14). The choice of method is 
dependent on experimental conditions 
including sample geometry and whether 
additional Investigation is required 
with the sample. For example, the 
dry-weight dimensional method can only 
be used for regularly-shaped samples 
with uniform dimensions such as cubes 
and rods, while liquid Immersion is 
potentially destructive due to liquid 
infusion into the sample. Because of 
the lack of a truly universal porosity 
fraction measurement method, it seems 
worthwhile to consider additional 
measurement/estimation methods that 
may be useful and convenient in cer- 
tain laboratory and industrial situa- 
tions. In this study, we consider the 
ultrasonic velocity measurement method 
for estimating porosity fraction. 

Ultrasonic velocity is a rela- 
tively simple measurement that requires 
the material specimen to have one pair 
of sides flat and parallel (ref. 15). 

The advantages of this method are that 
it is nondestructive and measurements 
can be made on different regions of a 
single specimen. Smith (ref. 2) and 
Nagarajan (ref. 16) were two of the 
first researchers to establish empiri- 
cal correlations between porosity frac- 
tion and ultrasonic velocity for 
polycrystalline materials. The corre- 
lations appeared relatively linear over 


the porosity fraction ranges investi- 
gated. Smith's work concerned metallic 
samples while Nagarajan' s work concerned 
ceramic samples. Other researchers 
began to investigate similar correla- 
tions with different materials. Here, 
we review and statistically analyze 
these empirical correlations between 
ultrasonic velocity and porosity frac- 
tion for polycrystalline materials. 
Initially, a semi-empirical model is 
developed showing the origin of the 
linear relationship between ultrasonic 
velocity and porosity fraction. Then, 
scatter plots of velocity versus per- 
cent porosity data are shown for A 1 2 O 3 , 
CuO, MgO, porcelain-based ceramics, 

PZT, SiC, SI 3 N 4 , steel, tungsten, UO 2 , 
(Uo. 3()P U 0. 70^C, and YBa 2 Cu 307 _ x . Lin- 
ear regression analysis produced pre- 
dicted slope, intercept, correlation 
coefficient, level of significance, and 
confidence interval statistics for the 
data. Additionally, velocity values 
predicted from regression for fully- 
dense materials are compared with those 
calculated from elastic properties. 


SEMI-EMPIRICAL MODEL 

When there are no boundary effects 
present, the velocity of a longitudinal 
ultrasonic wave traveling in a solid is 
related to the elastic properties and 
density of the solid by (ref. 17): 

V = {[E(l - v)]/[p(l + v)(l - 2v)]} 1/2 

(1 ) 

where V, E, p, and v are the veloc- 
ity, elastic modulus, bulk density, and 
Poisson's ratio, respectively, of the 
material. (The velocity of a shear 
ultrasonic wave traveling in a solid is 
related to the elastic properties and 
density of the solid by: 

V = (E/[2 P (1 + v)]}l/2 (1(a))) 

An "apparent" modulus (ref. 3) for 
porous materials can be considered 
which depends on the porosity fraction. 
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(5) 


Several early empirical Investigations 
provided evidence that the modulus 
increases exponentially with decreasing 
porosity fraction according to 
(refs. 18 and 19): 

E = E 0 exp(-bP' ) (2) 

where E 0 Is the elastic modulus of a 
fully-dense (nonporous) material, b 
Is an empirically-determined constant 
related to pore shape, pore distribu- 
tion, and the ( rat1o of open-to-closed 
pores, and P* Is the porosity frac- 
tion. The use of equation (2) to eval- 
uate E 0 by extrapolation from fitted 
experimental data has sometimes resulted 
in large discrepancies between the 
extrapolated and observed values 
(ref. 20). An alternative to equation 
(2) has been suggested to describe the 
relationship between elastic modulus 
and porosity fraction (ref. 21): 

E = E 0 ( 1 - P ' )2n+l (3 ) 

where n, like b. Is an empirically- 
determined constant that depends on 
pore distribution and pore geometry 
factors . 1 

Porosity fraction, P ' , can be 
expressed as: 

P' = (1 - (p/p 0 >) (4) 

where p 0 Is the theoretical (nonpo- 
rous material) density. Rearranging 
equation (4) allows us to express bulk 
density as a function of porosity frac- 
tion: 


'Concerning the relationship 
between Poisson's ratio and porosity 
fraction, most of the limited studies 
of Poisson's ratio show it decreasing 
with increasing porosity fraction less 
rapidly than for elastic modulus 
(ref. 3). In this development, it is 
assumed that Poisson's ratio is inde- 
pendent of porosity fraction. 


P = P 0 <1 - P'> 

Substituting equations (3) and (5) 
into equation (1) allows velocity to be 
expressed as: 

V = V 0 (l - P' ) n (6) 

where V 0 is a constant for a given 
material equal to: 

V 0 = { [ E 0 ( 1 - v)]/[p 0 (l + v)(l - 2v)]} 1/2 

(7) 

V 0 is the velocity in a fully-dense 
(nonporous) material, i.e., the "theo- 
retical" velocity. (For shear waves: 

V 0 = {E 0 /C2p 0 (l + v)]}l/2 (7(a))) 

The general case for all n can be 
shown by expanding the right-hand side 
of equation (6) using the binomial 
theorem (ref. 22) so that: 

V = V 0 {1 + [n(-P')3 + [n(n - 1) 
x (-P')2/2!] +. . .+ Cn(n - 1) 
x . . .(n - k + 1)(-P ) k /k!]+. . .} 

( 8 ) 

From the ratio test, equation (8) is 
absolutely convergent for |P| < 1. 

Setting n = 1 in equations (3) 
and (8) results in good agreement for 
a number of materials over a wide 
porosity fraction range (0.1 < P 
< 0.7) (ref. 21). In this case, the 
right-hand side of equation (8) is 
reduced such that: 

V = V 0 ( 1 - P ' ) (9) 

Equation (9) shows a linear relation- 
ship between velocity and porosity 
fraction and is the basis for select- 
ing linear regression to analyze the 
empirical correlations reported in 
this study. 
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It Is sometimes convenient to dis- 
cuss the relationship between velocity 
and percent porosity, %P, where: 

IP = (P')IOO (10) 

Solving equation (10) for P' and sub- 
stituting into equation (9) gives: 

V = m(7.P) + V 0 (11) 

where 

m = -V o /100 (12) 

Equation (11) shows a linear rela- 
tionship between V and %P where m 
and V 0 are the slope and intercept, 
respectively. 

We can also define a "percent 
theoretical velocity," %TV, where: 

%TV = (V/V 0 ) 1 00 (13) 

Solving equation (13) for V and 
substituting into equation (11) gives: 

%TV = m'(%P) + 100 (14) 

where 

m' = (m) 100/V o (15) 

Equation (14) shows a linear rela- 
tionship between UM and %P where 
m' and 100 are the slope and inter- 
cept, respectively. Presenting the 
velocity versus porosity fraction rela- 
tionship in terms of equation (14) is 
essentially a normalization procedure 
in that the theoretical velocity of a 
material and the type of wave (longitu- 
dinal or shear) used in the velocity 
measurement are "removed" as variables. 
From the derivative of equation (14), 
the following quantity can be defined: 

(1/m‘) = ( A%P / A%TV ) (16) 

where A is "change in." 


DATA EXTRACTION PROCEDURE 

Almost all of the data presented 
in this study were obtained from 
previously-published studies. The 
policy employed was that all of the 
available data should be tabulated and 
analyzed. In most cases, the reference 
provided V versus p data, either in 
the form of a table or plot. In some 
cases, the reference provided percent 
theoretical density (7op 0 ) or P' 
values instead of p values. Where 
necessary, p and 7.p 0 values were 
converted to 7.P values with the aid 
of equations (4) and (10). 2 


STATISTICAL ANALYSIS METHOD 

Linear regression analysis and its 
associated statistics utilized in this 
study are briefly described in the next 
several paragraphs. The authors felt 
that a comprehensive set of statistics 
was necessary for this analysis because 
of the variation among data sets in the 
number of velocity measurements and the 
porosity fraction range over which 
those measurements were made. 

Linear regression analysis is con- 
cerned with the problem of predicting 
or estimating the value of a (depend- 
ent) variable (V and °/=T V in equa- 
tions (11) and (14), respectively) on 
the basis of another (independent) var- 
iable (% P in equations (11) and (14)). 
For the sake of simplicity, we have 
applied the classical regression model 
(ref. 23) which involves the following 
assumptions. V (and %TV) has been 
assumed to contain all the error while 
IP has been assumed to contain no 


2 In most cases, the references 
provided V and p data to 3 or 4 
significant figures. For the sake of 
uniformity, all data and subsequent 
calculations including statistical 
values are presented in this report to 
at most 3 significant figures. 
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error. 3 The variance in V (and %TV) 
has been assumed to be constant for all 
values of %P, and the distribution 
about V has been assumed to be normal 
with mean values lying exactly on the 
regression line. It also has been 
assumed that only one V value was 
measured at a particular %P. 

Linear regression analysis results 
in predicted slope (m and m' in 
eqs. (11) and (14), respectively) and 
intercept (V 0 in eq. (ID) values that 
describe the relationship between V 
(and %TV) and %P. The Pearson 
product-moment correlation coefficient 
and level of significance statistics 
describe the quality of the regression. 
The correlation coefficient measures 
the strength of the linear relationship 
for the sample data. The level of sig- 
nificance, determined by the number of 
data points and the value of the corre- 
lation coefficient, determines an 
acceptance or confidence region for the 
regression. A level of significance of 
0.025 corresponds to a 95 percent con- 
fidence region. The smaller (better) 
the level of significance, the lower 
the probability that the value of the 
correlation coefficient can be attrib- 
uted to chance. 

Confidence intervals for the pre- 
dicted slope, intercept, and mean 
velocity values (the mean of further 
velocity measurements obtained at some 


^An analysis assuming errors in 
both variables is significantly more 
complicated. For some data sets, the 
uncertainty in 7.P may in fact be com- 
parable to that of V. The total 
uncertainties in each of XP and V 
including experimental uncertainties, 
uncertainties in extracting data from 
plots, and different assumed values of 
p 0 are estimated to be less than 
5 percent in all cases. 


XP value) are also presented. The 
95 percent symmetric confidence inter- 
val was chosen for the analysis. 4 In 
practical terms, the 95 percent confi- 
dence interval means that in 95 percent 
of the cases, the true value of the 
parameter will fall within the calcula- 
ted interval . 


REVIEW AND ANALYSIS 

The review and statistical analy- 
sis are presented in Table I and fig- 
ures 1 to 61 . Most of the figures show 
scatter plots of V versus XP data 
for A1 2O3 (refs. 16, 24 to 27), CuO 
(ref. 28) MgO (ref. 29), porcelain- 
based ceramics (refs. 7, 30 and 31), 

PZT (ref. 8), SiC (refs. 26, 32 to 35), 
SI3N4 (refs. 36 to 38), steel (ref. 39), 
tungsten (ref. 2), UO2 (ref. 40), (U0.30 
Pu 0.7(pC (ref. 41), and YBa2Cu307_ x 
(refs. 11, 42 to 45). The table presents 
the linear regression statistics corre- 
sponding to the scatter plots. The 
95 percent confidence interval for the 
predicted slope and intercept values 
are presented in table I while the 
95 percent confidence interval for mean 
predicted velocity values is shown by 


4 The choice of a particular size 
confidence interval is "economic" rather 
than mathematical. It depends directly 
on the cost of an error, and hence on 
the frequency with which one can afford 
to be wrong. High confidence intervals 
lead to wide limits, and if these limits 
are too wide to be useful, the gap 
between them must be reduced either by 
accepting less confidence or by increas- 
ing the amount of data (ref. 23). 
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dashed lines on the scatter plots. ® 

The quantity (A%P/A%TV) Is provided for 
all plot lines In the corresponding 
table entries. (Note that this quan- 
tity also has a confidence Interval 
associated with It, the width of which 
Is similar to that for m'). 

Pertinent Information concerning 
the reference's study including ultra- 
sonic technique, measurement uncertain- 
ties,® ml crostructural anisotropy, 
material processing techniques, and 
velocity variation within specimens is 
also included in the table. A blank 


^Several issues concerning the 
95 percent confidence intervals for 
predicted intercept, slope, and mean 
velocity values need to be noted. 

First, the assumption of only one V 
value for a particular IP value Is a 
conservative assumption that we know is 
false for some of the data sets in this 
review (see "comments" in table I). 

This assumption tends to make the lim- 
its of the confidence interval wider 
(worse) than if the confidence interval 
was calculated based on the mean of 
several velocity measurement values at 
a particular 7.P value. Second, the 
confidence intervals for predicted 
slope and intercept may not appear 
exactly symmetric in table I due to 
the round off procedure. Third, the 
95 percent confidence interval for mean 
predicted velocity values is in most 
cases drawn (dashed line) over the 
entire %P range shown. In several 
cases, the interval extends beyond the 
%P range where velocity data exists. 

In these cases, the interval widens 
(worsens) as expected where no data 
exists . 

®In most cases, the experimental 
uncertainties in the velocity and den- 
sity measurements were provided by the 
reference. In the event that they were 
not, the uncertainties were estimated 
from the reference's description of 
samples and measurement techniques, and 
from our experience. 


table entry Indicates that the informa- 
tion was unavailable. The "comments" 
in the table give the number of data 
points for that particular reference 
and In some cases point out a major 
conclusion determined by the reference 
concerning the V versus IP data. 

The figures are organized as fol- 
lows. Figures 1 to 4, 6 to 9, 12 to 
15, 17 to 23, 25, 26, 28, 29, 33 to 36, 
38, 40 to 42, 44, 45, 48 to 50, 52 to 

54, 57 to 59, show scatter plots of V 
(and %TV) versus t P for a single 
reference's data. 7 Results of multiple 
investigations for a specific material 
were also combined and analyzed as one 
data set in figures 5, 10, 37, 43, 46, 

55, and 60. Additionally, plots com- 
paring predicted (in most cases) 
regression lines obtained for a spec- 
ific material from different investi- 
gations are given in figures 11, 16, 

24, 27, 30 to 32, 39, 47, 51 , 56, and 
61. Where applicable for a material, 
scatter plots of longitudinal wave 
velocity data are presented before 
plots of shear wave velocity data. For 
most scatter plots, the solid line 
drawn is the linear regression line 
determined from the least-squares tech- 
nique. For the plots with only two 
data points, a line is drawn through 
the points. In this case, the correla- 
tion coefficient, level of significance, 
and confidence interval statistics are 
not appl icable. 


7 Because the range of %P values 
for which velocity measurements were 
obtained varied from reference to ref- 
erence, the ranges shown on the hori- 
zontal and vertical axes of the plots 
differ from one to the next, i.e., the 
plots are not standardized. In some 
cases, the plot axes had to be adjusted 
to allow the presentation of the 
95 percent confidence interval for mean 
predicted velocity values. 
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DISCUSSION 

General Observations 

Correlation coefficients with mag- 
nitudes greater than 0.95 were obtained 
in 31 out of 42 cases. Levels of sig- 
nificance with magnitudes less than 0.025 
were obtained in 36 out of 42 cases. 

For longitudinal wave velocity, pre- 
dicted intercepts (V 0 ) ranged from 
0.443 cm/)js for unpoled PZT4 and 
unpoled PZT5 of reference 8 to 1.23 cm/ 
ps for SIC of reference 33. For shear 
wave velocity, predicted intercepts 
(V 0 ) ranged from 0.313 cm/ps for YBa 2 
Cu 3 C> 7 _x to 0.786 cm/ps for SiC of ref- 
erence 26. The quantity (A%P/A%TV) 
ranged from -0.52 for porcelain of ref- 
erence 7 and poled PZT4 of reference 8 
to -8.26 for porcelain T2 of refer- 
ence 31. It is understandable that 
these quantities vary from one material 
to the next since each material has 
different elastic properties and density 
(eq. (1)). Predicted intercepts (V 0 ) 
for a specific material from different 
investigations agree fairly well (see 
the plots and tabular entries for 
AI 2 O 3 , SiC, and S i 3 N 4 ) . Predicted 
slopes for a specific material from 
different investigations agree fairly 
well in most cases. However, signifi- 
cant slope disparity is evident for 
AI 2 O 3 ; this may be due to the limited 
porosity range for the data of refer- 
ences 26 and 27 and the inclusion of 
green and prefired sample data in the 
cases of references 24 and 25. For 
reference 26, the limited porosity 
range over which data was obtained is 
manifested in extremely wide 95 percent 
confidence limits for predicted slope, 
intercept and mean velocity values. In 
fact, one of the bounds for the confi- 
dence limits for predicted slope is a 
positive value. 

Table II compares V 0 predicted 
from regression analysis with that cal- 
culated from equations (7) (longitudi- 
nal wave velocity) and (7a) (shear wave 
velocity) for several materials. 

Values of elastic modulus, Poisson's 


ratio, and density for fully-dense 
(single crystal and/or polycrystall 1 ne) 
materials used in the calculation are 
presented. The values of V 0 pre- 
dicted from regression and those 
obtained from calculation agree within 
approximately 17 percent in 16 out of 
16 cases, and within approximately 
6 percent in 11 out of 16 cases. 


Other Mi crostructural Variables 
Affecting Velocity 

Although porosity fraction seems 
to be a significant and perhaps the 
major mi crostructural feature affecting 
ultrasonic velocity, several references 
point to other mi crostructural varia- 
bles having an impact on velocity. 

These include slight compositional var- 
iations (ref. 31), preferred domain 
orientation (ref. 8 ), particle contact 
anisotropy (ref. 39), pore size distri- 
bution and geometry (ref. 2 ), and type 
of agglomeration (ref. 25). These var- 
iables may result in differences in 
predicted intercept (V 0 ) and slope for 
what is believed to be the same mate- 
rial from different investigations. 
Thus, the authors feel that the most 
accurate and precise application of the 
ultrasonic velocity method for estimat- 
ing porosity fraction first requires 
the development of accurate velocity 
versus porosity fraction relationships/ 
calibrations for the specific material 
of interest. 


Ramifications 

The estimation of batch-to-batch , 
sampl e-to-sampl e and wi thi n-sampl e %P 
variations for a material can be accom- 
plished if the quantity (A%P/A%TV) is 
known with reasonable confidence for 
that material. The nondestructive 
mapping of spatial porosity fraction 
variations within a sample by means of 
an ultrasonic scanning technique has 
been reported recently (refs. 50 and 
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51). This approach may also be useful 
in the analysis of the uniformity of 
composite materials (ref. 34). 


CONCLUSION 

A review and statistical analysis 
of the ultrasonic velocity method for 
estimating the porosity fraction in 
polycrystalline materials is presented. 
First, a semi-empirical model was deve- 
loped showing the origin of the linear 
relationship between ultrasonic veloc- 
ity and porosity fraction. Then, from 
a compilation of data produced by many 
researchers, scatter plots of velocity 
versus percent porosity data were shown 
for AI 2 O 3 , CuO, MgO, porcelain-based 
ceramics, PZT, SiC, Si 3 N 4 , steel, tung- 
sten, UO 2 , (Uo.3C)P u 0.7(pC, an d 
YBa 2 Cu 307 _ x . Linear ^egression analy- 
sis produced slope, intercept, correla- 
tion coefficient, level of significance, 
and confidence interval statistics for 
the data. Velocity values predicted 
from regression analysis for fully- 
dense materials are in good agreement 
with those calculated from elastic 
properties. The estimation of batch- 
to-batch, sample-to-sample, and within- 
sample variations in porosity fraction 
for a material can be accomplished with 
ultrasonic velocity measurements if 
reasonable confidence exists in the 
velocity versus percent porosity linear 
relationship. 


REFERENCES 

1. Goetzel, C.: Treatise on Powder 

Metallurgy. Wiley, 1963. 

2. Smith, J.T.; and Lopilato, S.A.: 
The Correlation of Density of 
Porous Tungsten Billets and 
Ultrasonic-Wave Velocity. Trans. 
Metal 1 . Soc. AIME, vol . 236, 

no. 4, Apr. 1966, pp. 597-598. 


3. Rice, R.W.: Microstructure Depend- 

ence of Mechanical Behavior of 
Ceramics. Properties of Micro- 
structure: Treatise on Materials 

Science and Technology, vol. 11 , 
R.K. MacCrone, ed., Academic Press, 
1977, pp. 199-381. 

4. Knudsen, F.P.; Parker, H.S.; and 

Burdick, M.D.: Flexural Strength 

of Specimens Prepared from Several 
Uranium Dioxide Powders— Its 
Dependence on Porosity and Grain 
Size Anol the Influence of Addi- 
tions of Titania. J. Am. Ceram. 
Soc., vol. 43, no. 12, Dec. 1960, 
pp. 641-647. 

5. Burdick, M.D.; and Parker, H.S.: 
Effect of Particle Size Bulk Den- 
sity and Strength Properties of 
Vranium Dioxide Specimens. J. Am. 
Ceram. Soc., vol. 39, no. 5, May 
1956, pp. 181-187. 

6 . Williams, E.C.; Reld-Jones, R.C.; 
and Dorrill, D.T.: Trans. Brit. 
Ceram. Soc., vol. 62, 1963, 

pp. 405-420. 

7. Boisson, J.; Platon, F.; and 

Boch, P.: Constanti Electriche E. 

Attrito Interno in Alconi Mater- 
iali in Funzione Della Porosita. 
Ceramurgia, anno VI, no. 2, 1976, 
pp. 74-76. 

8 . Patel, N.D.; and Nicholson, P.S.: 
Comparison of Piezoelectric Proper- 
ties of Hot-Pressed and Sintered 
PZT. Am. Ceram. Soc. Bull., 

vol . 65, no. 5, May 1986, 
pp. 783-787. 

9. Alford, N. McN.; et al . : The 

Effect of Density of Critical Cur- 
rent and Oxygen Stoichiometry of 

Y Ba 2 Cuz 0( X ) Superconductors. 
Nature, vol. 332, Mar. 3, 1988, 
pp. 58-59. 


8 



10. Song, S.N.; et al . : Hlgh-Tc Super- 

conductivity in Y-Ba-Cu-0 System. 
Adv. Ceram. Mater., vol . 2, no. 3B, 
1987, pp. 480-486. 

11. Blendel 1 , J.E.; et al.: 
Processing-Property Relations for 
YBZ 2 CU 3 0( 7 _ x ) H i gh-Tc Superconduc- 
tors. Adv. Ceram. Mater., vol. 2, 
no. 3B, July 1987, pp. 512-529. 

12. Jones, J.T.; and Berard, M.F.: 
Ceramics Industrial Processing and 
Testing, Iowa State University 
Press, 1972. 

13. Vander Voort, G.F.: Metallography 

Principles and Practices. 
McGraw-Hill , 1984, p.426. 

14. Clark, G.L.; and Liu, C.H.: Quani- 

titave Determination of Porosity by 
X-Ray Absorption. Anal. Chem., 
vol. 29, no. 10, Oct. 1957, 

pp. 1539-1541. 

15. Papadakls, E.P.: Absolute Accuracy 

of the Pulse-Echo Overlap Method 
and the Pulse Superposition Method 
for Ultrasonic Velocity. J. Acous. 
Soc. Am., vol. 52, no. 3, Pt. 2, 
Sept. 1972, pp. 843-846. 

16. Nagarajan, A.: Ultrasonic Study of 

Elasticity-Porosity Relationship in 
Polycrystalline Alumina. J. Appl . 
Phys., vol. 42, no. 10, Sept. 1971, 
pp. 3693-3696. 

17. Szilard, J.: Physical Principles 

of Ultrasonic Testing. Ultrasonic 
Testing, J. Szilard, ed., Wiley, 
1982, pp. 1-23. 

18. Ryshkewitch, E.: Compressive 

Strength of Porous Sintered 
Alumina and Zirconia to Ceramo- 
graphy. J. Am. Ceram. Soc., 
vol. 36, no. 2, Feb. 1953, 

pp. 65—68 . 


19. Spriggs, R.M.: Expression for 

Effect of Porosity on Elastic Modu- 
lus of Polycrystalline Refractory 
Materials, Particularly Aluminum 
Oxide. J. Am. Ceram. Soc., 

vol . 44, no. 12, Dec. 1961 , 
pp. 628-629. 

20. Soroka, I.; and Sereda, P.J.: 
Interrelation of Hardness Modulus 
of Elasticity and Porosity in Vari- 
ous Gypsum Systems. J. Am. Ceram. 
Soc., vol. 51, no. 6 , June 1968, 
pp. 337-338. 

21. Phani , K.K.; and Niyogi, S.K.: 
Porosity Dependence of Ultrasonic 
Velocity and Elastic Modulus in 
Sintered Uranium Dioxide — A 
Discussion. J. Mater. Sci. Lett., 
vol. 5, Apr. 1986, pp. 427-430. 

22. Swokowski , E.W.: Calculus With 

Analytical Geometry. Prindle, 

Weber & Schmidt, 1 975, p. 461 . 

23. Acton, F.S.: Analysis of Straight 

Line Data. Wi ley, 1959. 

24. Claytor, T.N.; et al . : Nondestruc- 

tive Measurement of Microstructure 
Evolution in Ceramics. Mater. 

Eval., vol. 47, no. 5, May 1989, 
pp. 532-535. 

25. Jones, M.P.; Blessing, G.V.; and 

Robbins, C.R.: Dry-Coupled Ultra- 

sonic Elasticity Measurements of 
Sintered Ceramics. Mater. Eval., 
vol . 44, no. 1 , June 1 986 , 

pp. 859-862. 

26. Stang, D.S.: Unpublished research, 

NASA Lewis Research Center, Cleve- 
land, Ohio. 

27. Arons, R.M.; and Kupperman, D.S.: 
Use of Sound-Velocity Measurements 
to Evaluate the Effect of Hot Iso- 
static Pressing on the Porosity of 
Ceramic Solids. Mater. Eval., 
vol. 40, no. 10, Sept. 1982, 

pp. 1076-1078. 


9 



28. Roth, D.J.; and Dolhert, L.E.: 
Unpublished research, NASA Lewis 
Research Center, Cleveland, Ohio, 
1990. 

29. D.S. Kupperman; and H.B. Karplus: 
Ultrasonic Wave Propagation Char- 
acteristics of Green Ceramics. Am. 
Ceram. Soc. Bull., vol . 63, no. 12, 
Dec. (1984), pp. 1505-1509. 

30. Prof. Dr. inz. Leszek Fi 1 i pczynski ; 
Dr. inz Zozislaw Pawlowski ; and 
Dr. inz Jerry Wehr, Ultrasonic 
Methods of Testing Materials. 2nd 
Edition, Translation, Butterworths , 
1966, p. 255. 

31. Shyuller, K.G.; Khennike, G.V.; and 

Kovziridze, Z.D.: Specific Fea- 

tures of Structure Evaluation in a 
Low-Temperature Ceramics. Glass 
Ceramics, vol. 44, 1988, 

pp. 473-475. 

32. Friedman, W.D.; et al . : Character- 

ization of Green Ceramics With 
X-Ray Tomography and Ultrasonics. 
Nondestructive Testing: High Per- 
formance Ceramics; Proceedings. 
American Ceramic Society, 1987, 

pp. 128-131. 

33. Baaklini, G.Y.; et al . : High Fre- 

quency Ultrasonic Characterization 
of Sintered Silicon Carbide. Am. 
Ceram. Soc., vol. 72, no. 3, March 
1989, pp. 383-387. 

34. Gruber, J.J.; Smith, J.M.; and 

Brockelman, R.H.: Ultrasonic 

Velocity C-Scans for Ceramic and 
Composite Material Characteriza- 
tion. Mater. Eval., vol. 46, 
no. 1, Jan. 1988, pp. 90-96. 

35. Kllma, S.J.; et al . : Ultrasonic 

Velocity for Estimating Density of 
Structural Ceramics. NASA 
TM-82765 , 1981. 


36. Derkacs, T.; Matay, I.M.; 

Brentnall, W.D.: Nondestructive 

Evaluation of Ceramics. Contract 
Report No.: AD-A027357, 

TRW-ER-7798-F ; TRW, Inc., May 1976. 

37. McLean, A.F., et al.: Brittle 

Materials Design, High Temperature 
Gas Turbine, Interim Report Num- 
ber 7, July 1, 1974 to December 31, 
1974, AMMRC CTR 75-8, Apr. 1975, 

p. 97. NOTE: RESTRICTED MATERIAL 

38. Thorp, J.S.; and Bushel 1 , T.G.: 
Ultrasonic Examination of Reaction- 
Bonded Silicon-Nitride. J. Mater. 
Sci., vol. 20, no. 6, June 1985, 
pp. 2265-2274. 

39. Papadakis, E.P.; and Petersen, 

B.W.: Ultrasonic Velocity as a 

Predictor of Density in Sintered 
Powder Metal Parts. Mater. Eval., 
vol . 37, no. 5, Apr. 1979, 

pp. 76-80. 

40. Panakkal , J.P.; and Ghosh, J.K.: 
Ultrasonic Velocity in Sintered 
Uranium-Dioxide Pellets. J. Mater. 
Sci. Lett., vol. 3, no. 9, 1984, 
pp. 835-836. 

41. Ghosh, J.K.; Panakkal, J.P.; and 

Roy, P.R.: Ultrasonic Velocity 

Measurement in The Characterisa- 
tion of Nuclear Fuel Pellets. 

World Conference on Nondestructive 
Testing, Butterworths, 1985, 

pp. 1201-1205. 

42. Gaiduk, A.L.; et al . : High- 

Frequency Acoustic Properties of 
Yttrium Ceramic. Sov. J. Low 
Temp. Phys., vol. 14, no. 7, July 
1988, pp. 395-402. 

43. Roth, D.J.; Dolhert, L.E.; and 

Kerwin, D.P.: Unpublished 

research, NASA Lewis Research 
Center, Cleveland, Ohio, 1990. 


10 



44. Ledbetter, H.M.; et al . : Elastic 

Constants and De BYe Temperature 
of Polycrystalline YBa 2 Cu 307 _ x - 
J. Mater. Res., vol . 2, no. 6, 
Nov/Dec. 1987, pp. 786-789. 

45. Round, R.; and Bridge, B.: 
Elastic-Constants of the High- 
Temperature Ceramic Superconductor 
YBa 2 Cu 307 _ x . J. Mater. Sci. Lett., 
vol . 6, no. 12, Dec. 1987, 

pp. 1471-1472. 

46. Alford, N.McN., et al.: Physical 

and Mechanical Properties of 
YBa 2 Cu 30 ( 7 _si > Superconductors. 
Mater. Sci., vol. 23, no. 3, March 
1988, pp. 761-768. 

47. G.E. Deiter, Mechanical Metallurgy. 
Third ed., McGraw-Hill, 1986, p.50. 

48. Engineering Property Data On 
Selected Ceramics. Vol. Ill, Single 
Oxides. Metals and Ceramics Infor- 
mation Center, Battel le's Columbus 
Laboratories, Columbus, OH, 1981, 
pp. 5.4.1-25. 


49. Engineering Property Data On 
Selected Ceramics. Vol. II., Car- 
bides. Metals and Ceramics Infor- 
mation Center, Battel 1 e * Columbus 
Laboratories, Columbus, OH, 1979, 
pp. 5.2.3-10. 

50. Generazio, E.R.; Roth, D.J.; and 

Baaklini, G.Y.: Acoustic Imaging 

of Subtle Porosity Variations in 
Ceramics. Mater. Eval . , vol. 46, 
no. 9, Sept. 1988, pp. 1338-1343. 

51. E.R. Generazio; D.J. Roth; and 

D.B. Stang: Ultrasonic Imaging of 

Porosity Variations Produced Dur- 
ing Sintering. J. Am. Ceram. 

Soc., vol. 72, no. 7, July 1989, 
pp. 1282-1289. 


11 



LIST OF SYMBOLS AND ABBREVIATIONS FOR TABLE I. 


V = Velocity (cm/^sec) 

SV = Shear wave velocity (cm//xsec) 

LV = Longitudinal wave velocity (cm/^sec) 

%TV = Percent theoretical velocity 
%P = Percent porosity 
Wt. = Weight 

( 3 q = Predicted value of intercept (Theoretical velocity) 

0' o = Predicted value of intercept (Percent theoretical velocity) 

= Predicted value of slope (Velocity/percent porosity) 

= Predicted value of slope (Percent theoretical velocity/percent porosity) 
N/A = Not applicable 
RBSN = Reaction-bonded silicon nitride 
A = Change in 

(A blank appearing in a table entry indicates that the author did not mention 
the subject or the information was otherwise unavailable.) 



TABLE I— ULTRASONIC VELOCITY VERSUS P 


Referent* Processing Nolet/ MIc restructure! 

Chemical Additives Anisotropy 


Velocity Velocity Density Density 

Measurement Measurement Measurement Measurement 

Technique Uncertainty Technique Uncertainty 

m (*) 


Theoretical 
Density, 
Used lo 
Calculate 
H Porosity 
(g/cc| 




Putse-echo 

overlap, 

longitudinal waves, 
2 2 5 MHz 


Dry-wt 

dimensional 


LV * -0 014 • %P + 0 817 


FOLDOUT FRAME / 

















































PERCENT POROSITY: REVIEW AND STATISTICAL ANALYSIS 



Predicted Une Equation 

Correlation 

95H Confidence 

95H Confidence 

AH P 

Urged 

Comments 

Corresponding 


(HTV * 0', »HP + 0 J 

Coefficient 

Intervals lor Predicted 

Interval! tor Predicted 

AHTV 

Velocity 


Figure Number 


level of 

Intercept (/ty and 

Intercept OS'!) and 


Variation 




Significance 

Slope OS,) 

Slope (0',) 


Across Ono 
Surface el 
Specimen 










<H) 



1 09 

HTV - -1 65 • %P + 100 

-0 999 

1 07 < 0 O £ 1 12 

97.7 £ 0' o £ 102 

-0 61 

15 

5 data points lumped together in low 

1 



0 0001 

-0019 £ 0, £ -0017 

-1.77 £ 0', £ -1.52 



H porosity region separated widely 
from 1 data point in high H porosrty 
{Partially-fired specimen) 


1 13 

HT V « -1.64 * HP + 100 

-0 992 

0 994 £ 0 Q £ 1 27 

87 8 £ 0' o s 112 

-0 81 

1 

Velocity may be very slightly sensitive 

2 



0 0001 

-0.022 < 0, s -0 015 

-1.94 s 3', fi - 134 



to the type of agglomerates found in 
ceramic samples, which may depend on 
whether sample is calcined or not 2 data 
points in low H porosity region widely 
separated from other 4 data points in 
high H porosity region (Green samples) 


1 00 

HTV - -0 736 • HP + 100 

-0 928 

0 987 £ 0 O £ 1 02 

98 3 £ 3'o * 102 

-1 36 


16 data points 

3 



0 0001 

-0 008 £ 0, s -0 007 

-0818 £ 0', £ -0 654 





1 01 

HTV - -0 377 • HP + 100 

-0.698 

0 923 £ 0 O £ 1 09 

92 2 < 0' o £ 108 

-2 65 

1 

Limited data region, 4 data points 

4 



0 302 

-0 016 < 0, < 0 008 

-1.54 £ 0’, ^ 0 782 





1 10 

HTV . - 1.43 • HP + 100 

-0 949 

1.06 s 3„ s 1 15 

95.7 £ 0' o £ 104 

-0 699 


32 data points, all longitudinal wave data 

5 



0.0001 

-0018 £ 3, £ -0013 

-1.61 £ 0’, £ -125 





) 669 

HTV « -0 477 • HP + 100 

-094 

0 646 £ 0 Q s 0 692 

96 7 £ (8' 0 < 103 

-2.10 

2 

6 data points 

6 



0 006 

-0 005 £ 0, £ - 0 002 

-0 669 £ 0', £ -0 284 





1655 

HTV - -1.87 • HP 4- 100 

-to 

0 648 < 0 Q < 0 663 

98.9 £ 0' o < 101 

-0 53 

15 

5 data points lumped together in low 

7 



0 0001 





H porosity region separated widely 





-0013 < 0, < -0012 

-1 93 £ 0', < - 1 81 



from 1 data point in high % porosity 
region (Partially-ftred specimen) 


3666 

HTV « -1 55 • HP + 100 

-0 987 

0 571 £ 0 O £ 0 762 

85 7 £ 0' o £ 114 

-0 65 

1 

Velocity may be very slightly sensitive 

8 



0 0003 

-0 008 5 0, < -0 013 

-1.90 £ 0', s -1 20 



to the type of agglomerates found in 
ceramic samples, which may depend on 
whether sample ts calcined or not, 2 data 
points in low H porosity region wxlely 
separated from other 4 data points m 
high H porosrty region (Green samples) 


628 

HTV - -0 662 • HP + 100 

-0 990 

0 622 < 0 O £ 0 635 

99 0 £ 0 q £ 101 

- 1 51 


17 data points 

9 



0 0001 

-0 005 £ 0, < -0 004 

-0 705 £ 0 ' | £ -0.610 





693 

HTV - -1 35 • HP + 100 

-0 910 

0 656 < d 0 < 0 729 

94 6 £ 0' o < 105 

-0.741 


35 data points, all shear wave data 

10 



0 0001 

-0011 £ 0, £ -0 006 

-157 £ 0', £ -1.14 












Comparison 

11 

474 

HTV » - 1 34 • HP + 100 

-0990 

0 400 < 0 Q < 0 549 

84 3 < 0' Q < 116 

-0.746 


4 data points, 4 specimens cut from 

12 



0 0100 





2 disks, CuO phase confirmed from 





- 0 009 < 0, < - 0 004 

-1 92 £ 0', < -0 760 



x-ray diffraction 


817 

N/A 

-1 0 

0673 s 0 O < 0961 

N/A 

-060 

2 

3 data points, matenal is unsintered. 

13 



0 012 





green compact 





-0 017 < 0, < -0 010 







FQLDQUI FRAME ^ 





















































































TABU I.— ULTRASONIC VELOCITY VERSUS PERCENT POROSITY 


Material 

Reference 

Processing Notes 
Chemical Additive 

J Microstructural 

* Anisotropy 

Velocity 

Measurement 

Technique 

Velocity 

Measuremen 

Uncertainty 

(*) 

Density 

t Measuremen 
Technique 

Density 

Measurement 

Uncertainty 

Theoretical 
Density. 
Used lo 
Calculate 
Ai Porosity 
<B/ec) 

Predicted 
line Equation 
(V = fi % • AVP + 

(A 

■ 

7 


Preferred granular 
orientation 

Pulse-echo 

overlap, 

Longitudinal waves 
2 MH ? 

<; 0 1 

Ory-wt 

dimensional 

s 1 

26 

IV - -0 013 • %P + 0 728 

%T 

H 

!■ 


Preferred granular 
orientation 

Pulse-echo 
overlap. Shear 
waves. 2 MH? 

s 0 1 

Dry-wl 

dimensional 

< 1 

■ 

■E 

I 


7 







■ 



B 

mm 



longitudinal waves, 
15-3 MH? 




■ 

LV - - 0 006 • %p + 0618 

%TV 

Porcelain Pi 
(See ref ) 

31 

(See ref.) 



B 



■ 


%TV 

Porcelain P2 
(See ref ) 

31 

(See ref ) 



< i 

'1 

B 

256 


%TV 

Porcelain T1 
(See ref.) 

31 

(See ref ) 



^ i 


B 

2 58 



IV - -0 001 • %P + 0 623 

%TV 

Porcelain T2 
(See ref.) 

31 

(See ref ) 


Longitudinal waves, 
15 MH? 

B 

1 


264 

LV * - 0 0008 *%P + 0 626 

%TV 

Unpoled PZT4 

6 


Insignificant 

■ 

■ 

— 

■ 

1 

80 

IV « - 0 007 • %P + 0 443 

At TV 

Poled PZT4 

8 


Preferred electrical 
domain orientation 

Pulse -echo 
overlap. 

Longitudinal waves. 
50 MH? 

■ 

H 


80 

IV « - 0 009 • %P + 0 483 


PZT4 

8 










Unpoled PZT5 

8 


Insignificant 

Pulse-echo 

overlap. 

Longitudinal waves. 
50 MH? 

■ 


■ 

80 

IV - -0 008 ♦ %P + 0 443 


Poled PZT5 

8 


Preferred electrical 
domain orientation 

Pulse-echo 

overlap, 

Longitudinal waves, 
50 MH? 

■ 


■ 

80 

IV = -0 010 • %P + 0 486 

%TV = 

PZT5 

8 











FOLDOUT FRAME J 













































POBOSJTY: REVIEW AND STATISTICAL ANALYSIS- CONTINUED 



Predicted Line Equation 
(*TV*0 W + 

Correlation 

Coefficient 

964fa Confidence 
Intervals for Predicted 
Intercept and 

Slope (0,) 

95At Confidence 
Intervals for Predicted 
Intercept (0' a ) and 
Slope (0 ,) 

Att P 
A* TV 

Largest 
Velocity 
Variation 
Across One 
Surlace of 
Specimen 

m 

Comments 

1 

Level of 
Significance 

728 

%TV - -1.78 • %P + 100 

-0 993 
0 0001 

0 713 5 0 O 5 0.743 
-0 0014 s ^ < -00012 

97 9 < 0' o * 102 
-1 87 s 0', < -1 70 

-056 


27 data pants, multiphase ceranvc, 
sample composition and thus theoretical 
density vary slightly with firing temper- 
ature; measurements made along 
extrusion axis 

14 

440 

%TV - -1 93 • %P ♦ 100 

-0998 
0 0001 

0.442 < 0 O £ 0 454 
-0 009 < 0, £ -0 008 

mum 

a 

o 

t 


15 








Comparison 

16 

618 

%TV - -0 935 • %P + 100 


0 606 < 0 O < 0 630 
-0,007 sj3, £ -0 005 

98 1 £ 0' o ^ 102 
-1.11 5 0', £ -0760 

-1.07 


7 data points, multiphase ceramic, 
sample composition and thus theoretical 
density vary slightly with firing temp 

17 

611 

<*TV - -0 312 • %P + 100 

- 0 586 
0 414 

0 553 ^ 0 Q < 0 668 
-0010 < 0, £ 0.006 

90 5 £ 0' o < 109 
-1 34 £ 0', < -1 ?8 

-321 


4 data points, multiphase ceramic, 
sample composition and thus theoretical 
density vary slightly with firing temp 

18 

615 

%TV - - 0 740 ♦ %P + 100 

-0 903 
0.017 

0.598 5 0 O £ 0.631 
-0 007 5 0, £ -0 002 

97 3 £ 0' o £ 102 
-1.04 ^ 0', £ -0 436 

-1 35 


19 

623 

i*) TV - -0 167 • P + 10Q 

-0 894 
0106 

0 612 < 0 O < 0 632 
-0 003 < 0, < 0 0005 

98 6 < 0' o £ 101 
-0.422 < 0', £ 0 088 

-5 99 


20 

626 

%TV - -0 121 • %P + 100 

-0 947 
0 053 

0 621 < 0 O < 0 631 
-0 002 < 0, < 0 003 

99 2 < 0' o £ 101 
-0 246 s 0', £ 0 004 

-8 26 


21 

443 

%TV - - 1 60 • %P + 100 

N/A 

N/A 

N/A 

-060 


2 data points for each set, PZT4, PZT5, 
and PZT7 are each different solid solution 
combinations of PbZr0 3 and PbTt0 3 but 
theoretical density is - 8 0 Q/cc for each, 
poling orients electric domains, velocity 
measured along polanzed direction for 

22 

N/A 

23 

403 

%TV - - 1.93 • %P ♦ 100 

N/A 

N/A 

N/A 

-0.52 


poled specimens; no correlation coef- 
ficient is given since only 2 data points 
for each set, ail poled sets have higher 
theoretical velocity than corresponding 
unpoled set 

N/A 

__ 







Companson 

24 

443 

*t>TV - - 1 85 • %P + 100 

N/A 

N/A 

N/A 

-0 54 


2 data points for each set. PZT4 PZT5, 
and PZT7 are each different solid solution 
combinations of PbZr0 3 and PbTt0 3 but 
theoretical density is -8.0 g/cc lor each, 
poling orients electric domains, velocity 
■ measured along polarized direction for 

25 

N/A 

26 

406 

%TV - -2 12 • %P + 100 

N/A 

N/A 

N/A 

- 0 47 


poled specimens, no correlation coef- 
ficient is given since only 2 data pants 
for each set, all poled sets have higher 
theoretical velocity than corresponding 
un poled set 

N/A 








Companson 

27 


fgldout FRAME 3, 


15 




































































26. 33-35 


36 



3 22 SV » - 0 009 • %P + 0 7 


Hot-pressed silicon 
nitride has 1% 
MgO sintering aid 
and 0 5-1% impu- 
tes, R8SN has 
< 1% impurities 
and various 
amounts of unre- 
acted silicon 


Insignificant 


I Pulse-echo 
overlap. 

longttudinal waves, 
25-45 MHz 


Ory-wt 

dimensional 


-COUT FRAME I 






















iT POROSITY: REVIEW AND STATISTICAL ANALYSIS-CONTINUED 



Predicted Une Equation 

Correlation 

95A* Confidence 

9544 Confidence 


Lirgest 

Comments 

Corresponding 


{%TV * ^ ' 1 • W + 0 J 

Coefficient 

Intervals for Predicted 

Intervals for Predicted 

TV 

Velocity 


Figure Number 

1 

level el 

Intercept Oty and 

Intercept (3 # ) end 


Variation 




Significance 

Slope I3,f 

Slope (3 ,1 


Across One 
Surlace ef 
Specimen 









m 



) 464 

%TV - -0 760 • %P + 100 

N/A 

N/A 

N/A 

-1 32 


2 data points for each set, PZT4. PZT5, 

28 



N/A 





and PZT7 are each different soIkJ solution 
combinations ot PbZr0 3 and Pb Tr0 3 but 
theoretical density is -8.0 g/cc for each; 
poling orients eleclnc domains, velodty 
measured along polanzed direction for 










0 494 

%TV - -0 952 • %P + 100 

N/A 

N/A 

N/A 

-1 05 


poled specimens, no correlation coef- 

29 








ficient is given since only 2 data points 
for each set, all poled sets have higher 
theoretical velocity than corresponding 
unpoied set 









Comparison 

30 








Companson 

31 

9 







Comparison 

32 

Q 

N/A 

- 0 974 

0 397 < 3 0 < 0 535 

N/A 

-0 633 

30 

8 data points, matenal is unsintered 

33 

I 


0 0001 

-0 009 s 3, < -0 006 




green compact, large velocity vanation 
even though small density vanation within 
specimen; quality of contact between 
individual powder particles may affect 
velocity 


23 

%TV » -1 16 • %P + 100 


1 22 < )3 0 < 1 24 


-0 862 

< 0.1 

8 data points, each data point is lor a 

34 



0.0001 

-0016 £ ^ 5 -0013 




particular batch and is the average of 
- 24 measurements on - 8 specimens 
Velocity not greatly sensitive to mean 
pore size, mean pore orientation and 
mean grain size 


21 

%TV * - 0 916 * %P ♦ 100 

-0 999 

1 21 s 3 0 s 1 22 


- 1 09 

< 5 

6 data points 

35 



0.0001 

-0012 < 3, < -0011 






22 

%TV « -0 883 • %P + 100 

-0 957 

1 22 £ 3 0 £ 1 22 


- 1 13 

< 1 

194 data points 

36 



0 0001 

-0011 s 3, < -0010 






22 

%TV - -0 912 • %P + 100 

-0 964 

1 216 s 3 0 < 1 220 


- 1 10 


208 data points, all longitudinal wave data 

37 



0 0001 

-001? < 3, s -0011 






786 

%TV - - 1.16 • %P + 100 

-0 991 

0 711 < 3 0 s 0 859 

990 s 3' 0 £ 101 

-086 

< 0 1 

8 data points, most specimens Irom 

38 



0 0001 

-0 010 < 3, S -0 008 

- 1 32 < 3', < -1 01 



same batches used by Baaklmi 









Comparison 

39 

11 

%TV - - 1.34 • %P «■ 100 

N/A 

N/A 

N/A 

-0 75 


No correlation coefficient is given since 

40 



N/A 





only 2 data points, but each point is 
average of measurement on 25 specimens, 
data point in low % porosity region is for 









hot-pressed S 3 N 4 while data point in high 
% porosity is for reaction- bonded $i 3 N 4 



fOLDOUf FRAME A 


17 




































































TABU I.— ULTRASONIC VELOC TTT VERSUS PERCENT : 



Material 

Reference 

Processing Notes/ 
Chemical Additives 

Mterostructural 

Anisotropy 

Velocity 

Measurement 

Technique 

Velocity 

Measurement 

Uncertainty 

m 

Density 

Measurement 

Technique 

Density 

Measurement 

Uncertainty 

m 

Theoretical 
Density, 
Used to 
Calculate 
H Porosity 
(g/cc) 

Predicted 
Line Equation 
(V * * AM* + tfj) 

s» 3 n 4 

37 

Injection- molded, 
slip-cast, and hot 
pressed specimens 


Thru-transmission 
transit time. 
Longitudinal waves, 
5 MHz 




3 30 

LV - -0 016 * HP + 1 ‘ 

St 3 N 4 

38 

RBSN has < 1% 
impurities and var- 
ious amounts of 
unreacted silicon 

Insignificant 

Pulse-echo 

overlap, 

Longitudinal waves, 
15 MHz 

s 0 1 

Dry-wt 

dimensional 

s 0 1 

330 

LV - -0 013 • %P + 1 - 

Si 3 N 4 

36-38 



Longitudinal waves 




3 30 

LV - -0.014 • HP ♦ 1 * 

Si 3 N 4 

37 

Injection-molded, 
slip-cast and hot 
pressed specimens 


Thru- tra ns mission 
transit time, Shear 
waves. 5 MHz 




3 30 

SV = -0.008 • HP ♦ 0 6: 

s. 3 n 4 

38 

RBSN has < 1% 
impurities and var- 
ious amounts of 
unreacted silicon 

Insignificant 

Pulse-echo 
overlap, Shear 
waves, 15 MHz 

< 0 1 

Dry-wt 

dimensional 

s 0 1 

330 

SV - -0 007 * %p + 0 6' 

^11 




Shear waves 




330 

SV = - 0 007 • HP ♦ 0 


s, 3 n 4 

36-38 










Steel 

A-direction 
{See ref ) 

39 


Possibly particle 
contact amstropy 
{ < 5%) based on 

Thru -transmission 
pulse-echo overlap. 
Dry coupling Longi- 
tudinal waves, 

1 5-2 25 MHz 

< 001 

ASTM 

B-328-60 

Undetermined 
probably 
< 2% 

7 85 

LV «= -0 007 • HP ♦ 0 5£ 


Steel 

B-direction 
(See ref } 

39 


pressing direction 

Thru-transmission 
pulse-echo overlap. 
Dry coupling Longi- 
tudinal waves. 

1 5-2 25 MHz 

< 0 01 

ASTM 

B-328-60 

Undetermined. 

probably 

7 85 

LV - - 0 009 • HP + 0 58 


B 




Thru-transmission 
pulse-echo overlap. 
Dry coupling Longi- 
tudinal waves. 

1 5-2 25 MHz 

< 001 

ASTM 

B-328-60 

Undetermined, 
probably 
< 2% 

7 85 

LV - -0 009 * HP + 0 59C 


Steel 

39 










FOLDOUj' rRAME, I 


























































a POROSmr: REVIEW and statistical analysis-continued 



I Predicted Line Equation 

Correlation 

95H Confidence 

95H Confidence 

aHP 

Largest 

Comments 



<HTV 

M , • + 

Coefficient 

Intervals lor Predicted 

Intervals for Predicted 

aH TV 

Velocity 





Level of 

Intercept (£,) and 

Intercept (0'g) and 


Variation 
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Significance 

Slope (0,) 

Slope (0 ' , ) 


Across One 
Surface of 
Specimen 
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13 data points 
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5 data points 
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dependence of velocity, velocity also may 
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TABLE L— ULTRASONIC VILOCITY VERSUS PERCENT POROSITY: RE 


Material 

Reference 

Processing Notes/ 
Chemical Additives 

Mlcroslructuril 

Anisotropy 

Velocity 

Measurement 

Technique 

Velocity 

Measurement 

Uncertainty 

m 

Density 

Measurement 

Technique 

Density 

Measurement 

Uncertainty 

m 

Theoretical 
Density, p t . 
Used to 
Calculate 
V* Porosity 

(fl/CC) 

Predicted 
Line Equation 

(V = V M i p + /»*) 

Predi 

(HTV 

Tungsten 

2 

Kenna nominal 
4 micron starling 
powder size 


Thru- transmission 
pulse-echo overlap. 
Longitudinal waves, 
1 MHz 



m 

193 

LV = - 0 005 • %P + 0 520 

%TV = 

Tungsten 

2 

G E. nominal 
4 micron starting 
powder size 


Thru-transmission 
pulse-echo overlap. 
Longitudinal waves. 
1 MHz 

s 0 5 

ASTM 

C20-46 

Liquul- 

immersion 


19 3 

LV - - 0 006 • %P + 0 558 

%TV = 

Tungsten 

2 

G E nominal 
18 micron starling 
powder size 


Thru-transmission 
pulse-echo overlap. 
Longitudinal waves. 
1 MHz 


ASTM 

C20-46 

Liquid 

immersion 


19 3 

LV = - 0 008 • %P + 0 554 

%TV = 

H 

2 



Longitudinal waves 




193 

LV = - 0 006 • %P + 0 533 

r 

%TV . 

Tungsten 

2 










UO ? 

40 



T hickness -cum- 
veiocity meter, 
Longitudinal waves, 
10 MHz 

■ 


■ 

10 96 

LV = - 0 008 • %P + 0 550 

%TV = 

W 0 30 Pu O 70^ 

41 



Thickness-cum- 
veloctty meter, 
Longitudinal waves. 
15 MHz 

< 1 


■ 

12 19 

LV = - 0 004 • %P 4 0 460 

%TV = - 

YBa 2 Cu 3 0 7 _ t 

43 

Different starting 
powders 


Puise-echo/cross- 
cor relation Longitu- 
dinal waves, Dry 
and wet coupling, 
5-20 MHz 

s 04 

Dry-wt 

dimensional 

■ 

638 

LV - - 0 007 • %P + 0 565 

%TV * - 
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POROSITY: REVIEW ANO STATISTICAL ANALYSIS-CONTINUED 


B 

Predicted Line Equation 

Correlation 

95H Confidence 

95H Confidence 

AH P 

Urges! 

Comments 

Corresponding 
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{HTV » 0 ' j • HP + 0’,) 

Coefficient 

Intervals for Predicted 
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Velocity 


Figure Number 


Level of 

Intercept 0$*) *nd 

Intercept (0 t ) and 


Variation 





Significance 

Slope (0,) 

Slope {0 ' , ) 


Across One 
Surface of 
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m 
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HTV - -0 939 • HP + 100 

-0 960 

0 465 < 0 O < 0.575 

89 4 5 0' fl s HI 

-1 06 


Velocity shown lo be sensitive to pore 

52 



0010 

-0 007 < 0 1 < -0 003 

-138 5 0', < -0 501 



$» 2 e distnbution/mean pore size and gram 
size distnbution/mean gram size which t$ 
a function of slartmg powder size disto- 
rtion. These results show (jiff erent theo- 
retical velocity for tungsten depending on 


558 

HTV - -1.13 • HP + 100 


0,387 s 0 Q < 0 729 

69 3 s 0' o s 131 

-088 


53 



0 028 





the starting powder size 





-0 007 ^ 3, £ -0 006 

-1.92 5 0', s -0 351 



5 data points for Kenna 4 powder. 

5 data points for GE 4 powder, 

4 data points for GE 18 powder 

Each data point is average of 17 meas- 
urements across sample 
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54 



0 008 
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14 data points, all longitudinal wave data 
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56 
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17 data points 
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4 data points 
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13 data points 
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TABLE I.— ULTRASONIC VELOCITY VERSUS PERCENT P 


Material 

Reference 

Processing Notes/ 
Chemical Additives 

Microstmctural 

Anisotropy 

Velocity 

Measurement 

Technique 

Velocity 

Measurement 

Uncertainty 
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Density 

Measurement 
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Measurement 
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Density, p,, 
Used to 
Calculate 
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(g/cc) 

Predicted 
Line Equation 
(V = • *P + /y 


43 
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Pulse-echo/cross- 
correlaticm, Shear 
waves, 10 MHz 

s 04 

Dry-wt 

dimensional 

< 1 


3B 

SV - - 0 002 

• %P + 0.31 

42 



Phase comparison 
method, Shear 
waves, 50 MHz 

< 2 



6 

YBa 2 Cu 3 0 7 _ t 
Superconductor 

11 



Pulse-echo overlap, 
Shear waves. 

3-10 MHz 

s 1 

dimensional 

Dry wt 

< 1 


44 



PJse-echo matchup. 
Shear waves, 

3-4 MHz 

H 

■ 

■ 


45 
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■ 

■ 

YB^CUjOy _ x 
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POROSITY: REVIEW AND STATISTICAL AHALYSIS-CONCLUOED 
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PERCENT POROSITY 

Figure 1 . - Longitudinal velocity versus percent porosity for Al 2 0 3 (ref. 24). 
Velocity - -0.018 x percent porosity + 1 .09. 

Percent theoretical velocity * -1.65 x percent porosity + 100. 

Correlation coefficient - -0.999. 
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PERCENT POROSITY 

Figure 2. - Longitudinal velocity versus percent porosity for Al 2 0 3 (ref. 25). 
Velocity * -0.019 x percent porosity + 1 .13. 

Percent theoretical velocity = -1 .64 x percent porosity + 100. 

Correlation coefficient = -0.992. 
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PERCENT POROSITY 

Figure 3. - Longitudinal velocity versus percent porosity for Al 2 0 3 (ref. 16). 
Velocity « -0.007 x percent porosity + 1 .00. 

Percent theoretical velocity - -0.736 x percent porosity + 100. 
Correlation coefficient - -0.982. 
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PERCENT POROSITY 


Figure 4. - Longitudinal velocity versus percent porosity for AI2O3 (ref. 26). 
Velocity = -0.004 x percent porosity + 1 .01 . 

Percent theoretical velocity = -0.377 x percent porosity + 100. 
Correlation coefficient = -0.698. 



PERCENT POROSITY 

Figure 5. - Longitudinal velocity versus percent porosity for AI2O3 (refs. 16, 24 to 26). 
Velocity * -0.016 x percent porosity + 1 .10. 

Percent theoretical velocity = -1 .43 x percent porosity + 100. 

Correlation coefficient = -0.949. 



PERCENT POROSITY 

Figure 6. - Shear velocity versus percent porosity for AI2O3 (ref. 27). 
Velocity = -0.003 x percent porosity + 0.669. 

Percent theoretical velocity = -0.477 x percent porosity + 100. 
Correlation coefficient = -0.936. 
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Figure 7. - Shear velocity versus percent porosity tor AI2O3 (ref. 24). 
Velocity « -0.01 2 x percent porosity + 0.655. 

Percent theoretical velocity « -1 .87 x percent porosity + 100. 
Correlation coefficient - -1 .00. 
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PERCENT POROSITY 


Figure 8. - Shear velocity versus percent porosity for AI2O3 (ref. 25). 
Velocity - -0.010 x percent porosity + 0.666. 

Percent theoretical velocity ■ -1 .55 x percent porosity + 100. 
Correlation coefficient - -0.987. 
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PERCENT POROSITY 


Figure 9. - Shear velocity versus percent porosity for AI2O3 (ref. 16). 
Velocity - -0.004 x percent porosity + 0.628. 

Percent theoretical velocity « -0.662 x percent porosity + 100. 
Correlation coefficient « -0.990. 
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Figure 10. - Shear velocity versus percent porosity for AI 2 O 3 (refs. 16, 24, 25, 27). 
Velocity -= -0.009 x percent porosity + 0.693. 

Percent theoretical velocity = -1 .35 x percent porosity + 100. 

Correlation coefficient = -0.910. 



Figure 1 1. - Ultrasonic velocity versus percent porosity for Al 2 0 3 (refs. 16, 24 to 27). 



Figure 12. - Longitudinal velocity versus percent porosity for CuO (ref. 28). 
Velocity = -0.006 x percent porosity + 0.474. 

Percent theoretical velocity = -1 .34 x percent porosity +100. 

Correlation coefficient = -0.990. 
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PERCENT POROSITY 

Figure 13. - Longitudinal velocity versus percent porosity for "green" M fl O (ref. 29). 
Velocity = -0.018 x percent porosity + 0.817. 

Correlation coefficient = -1 .00. 
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PERCENT POROSITY 

Figure 14. - Longitudinal velocity versus percent porosity for porcelain (ref. 7). 
Velocity = -0.013 x percent porosity + 0.728. 

Percent theoretical velocity = -1 .78 x percent porosity + 100. 

Correlation coefficient = -0.994. 



PERCENT POROSITY 

Figure 15. - Shear velocity versus percent porosity for porcelain (ref. 7). 
Velocity = -0.009 x percent porosity + 0.448. 

Percent theoretical velocity * -1 .93 x percent porosity + 100. 
Correlation coefficient - -0.998. 
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PERCENT POROSITY 


Figure 16. - Ultrasonic velocity versus percent porosity for porcelain (ref. 7). 
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PERCENT POROSITY 

Figure 17. - Longitudinal velocity versus percent porosity for porcelain (ref. 30). 
Velocity = -0.006 x percent porosity + 0.618. 

Percent theoretical velocity - -0.935 x percent porosity + 100. 

Correlation coefficient = -0.987. 



PERCENT POROSITY 

Figure 18. - Longitudinal velocity versus percent porosity for porcelain PI (ref. 31). 
Velocity = -0.002 x percent porosity + 0.61 1 . 

Percent theoretical velocity * -0.312 x percent porosity + 100. 

Correlation coefficient = -0.586. 
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PERCENT POROSITY 

Figure 19. * Longitudinal velocity versus percent porosity for porcelain P2 (ref. 31). 
Velocity - -0.005 x percent porosity + 0.61 5. 

Percent theoretical velocity - -0.740 x percent porosity + 100. 

Correlation coefficient = -0.983. 
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PERCENT POROSITY 

Figure 20. - Longitudinal velocity versus percent porosity for porcelain T1 (ref. 3). 
Velocity - -0.001 x percent porosity + 0.623. 

Percent theoretical velocity = -0.167 x percent porosity + 100. 

Correlation coefficient = -0.894. 



PERCENT POROSITY 

Figure 21 . * Longitudinal velocity versus percent porosity for porcelain T2 (ref. 31 ). 
Velocity = -0.001 x percent porosity + 0.626. 

Percent theoretical velocity = -0.121 x percent porosity + 100. 

Correlation coefficient = -0.947. 
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PERCENT POROSITY 


Figure 22. - Longitudinal velocity versus percent porosity for unpoled PZT4 (ref. 8). 
Velocity = -0.007 x percent porosity + 0.443. 

Percent theoretical velocity = -1 .68 x percent porosity + 100, 



Figure 23. - Longitudinal velocity versus percent porosity for poled PZT4 (ref. 8). 
Velocity = -0.009 x percent porosity + 0.483. 

Percent theoretical velocity = -1 .93 x percent porosity + 100. 
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Figure 24. * Ultrasonic velocity versus percent porosity for PZT4 (ref. 8). 
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PERCENT POROSITY 


Figure 25. - Longitudinal velocity versus percent porosity for unpoled PZT5 (ref. 8). 
Velocity « -0.008 x percent porosity + 0.443. 

Percent theoretical velocity * -1 .85 x percent porosity + 100. 



Figure 26. - Longitudinal velocity versus percent porosity for poled PZT5 (ref. 8). 
Velocity « -0.010 x percent porosity + 0.486. 

Percent theoretical velocity - -2.12 x percent porosity + 100. 



Figure 27. - Ultrasonic velocity versus percent porosity for PZT5 (ref. 8). 
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Figur© 28. - Longitudinal velocity versus percent porosity tor unpoled PZT7 (ref. 8). 
Velocity = -0.004 x percent porosity + 0.464. 

Percent theoretical velocity - -0.760 x percent porosity + 100. 
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Figure 29. - Longitudinal velocity versus percent porosity for poled PZT7 (ref. 8). 
Velocity * -0.005 x percent porosity + 0.494. 

Percent theoretical velocity * -0.952 x percent porosity + 100. 



Figure 30. - Ultrasonic velocity versus percent porosity for PZT7 (ref. 8). 
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Figure 31 . - Ultrasonic velocity versus percent porosity for unpoled PZT (ref. 8). 



Figure 32. - Ultrasonic velocity versus percent porosity for poled PZT {ref. 8). 



Figure 33. - Longitudinal velocity versus percent porosity for green SiC (ref. 32). 
Velocity - -0.007 x percent porosity + 0.464. 

Correlation coefficient = -0.974. 
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Figure 34. - Longitudinal velocity versus percent porosity for SiC (ref. 33). 
Velocity = -0 .0 1 4 x percent porosity + 1 .23. 

Percent theoretical velocity = -1 .16 x percent porosity + 100. 
Correlation coefficient = -0.993. 



Figure 35. - Longitudinal velocity versus percent porosity for SiC (ref. 34). 
Velocity = —0.01 1 x percent porosity + 1 .21 . 

Percent theoretical velocity = -0.916 x percent porosity + 100. 
Correlation coefficient = -0.999. 
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Figure 36. - Longitudinal velocity versus percent porosity for SiC (ref. 35). 
Velocity = -0.01 1 x percent porosity + 1 .22. 

Percent theoretical velocity - -0.883 x percent porosity + 100. 
Correlation coefficient = -0.957. 
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PERCENT POROSfTY 

Figure 37. - Longitudinal velocity versus percent porosity for SiC (refs. 33 to 35). 
Velocity - -0.01 1 x percent porosity + 1 .22. 

Percent theoretical velocity « -0.912 x percent porosity + 100. 

Correlation coefficient - -0.964. 



Figure 38. - Shear velocity versus percent porosity for SIC (ref. 26). 
Velocity * -0.009 x percent porosity + 0.786. 

Percent theoretical velocity » -1 .16 x percent porosity + 100. 

Cor relation coefficlen t = -0 . 99 1 . 



Figure 39. - Ultrasonic velocity versus percent porosity for SIC (refs. 26, 33 to 35). 
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Figure 40. - Longitudinal velocity versus percent porosity for Si 3 N 4 (ref. 36). 
Velocity = -0.01 5 x percent porosity +1.11. 

Percent theoretical velocity = -1.34 x percent porosity + 100. 
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Figure 41 . - Longitudinal velocity versus percent porosity for Si 3 N 4 (ref. 37). 
Velocity = -0.016 x percent porosity + 1.14. 

Percent theoretical velocity = -1 .41 x percent porosity + 100 
Correlation coefficient = -0.997. 



PERCENT POROSITY 

Figure 42. - Longitudinal velocity versus percent porosity for Si 3 N 4 (ref. 38). 
Velocity = -0.013 x percent porosity + 1.12. 

Percent theoretical velocity = -1 .1 5 x percent porosity + 100 
Correlation coefficient = -0.991 . 
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Figure 43. - Longitudinal velocity versus percent porosity for SI 3 N 4 (refs. 36 to 38). 
Velocity « -0.014 x percent porosity + 1 .12. 

Percent theoretical velocity * -1 .27 x percent porosity + 100 
Correlation coefflclen t * -0.981 . 



Figure 44. - Shear velocity versus percent porosity for (ref. 37). 
Velocity - -0.008 x percent porosity + 0.652. 

Percent theoretical velocity - -1 ,18 x percent porosity + 100 
Correlation coefficient * -0.991 . 
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Percent theoretical velocity - -1 .10 x percent porosity + 100 
Correlation coefficient « -0.984. 
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Figure 46. - Shear velocity versus percent porosity for Si 3 N 4 (refs. 37, 38). 
Velocity = -0.007 x percent porosity + 0.645. 

Percent theoretical velocity = -1 .07 x percent porosity + 100 
Correlation coefficient = -0.973. 



Figure 47. - Ultrasonic velocity versus percent porosity for Si 3 N 4 (refs. 36 to 38). 



Figure 48. - Longitudinal velocity measured in the A direction versus percent 
porosity for steel (ref. 39). 

Velocity = -0.007 x percent porosity + 0.563. 

Percent theoretical velocity =-1 .19 x percent porosity + 100. 

Correlation coefficient = -0.972. 
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Figure 49. - Longitudinal velocity measured in the B direction versus percent 
porosity for steel {ref. 39). 

Velocity - -0.009 x percent porosity + 0.588. 

Percent theoretical velocity - -1 .53 x percent porosity + 100. 

Correlation coefficient = -0.985. 



Figure 50. - Longitudinal velocity measured in the C direction versus percent 
porosity for steel (ref. 39). 

Velocity * -0.009 x percent porosity + 0.590. 

Percent theoretical velocity = -1 .48 x percent porosity + 100. 

Correlation coefficient = -0.996. 



Figure 51 . - Ultrasonic velocity versus percent porosity for steel (ref. 39). 
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Figure 52. - Longitudinal velocity versus percent porosity tor tungsten (ref. 2) with 
Kenna nominal 4 pm starting powder size. 

Velocity = -0.005 x percent porosity + 0.520. 

Percent theoretical velocity = -0.939 x percent porosity + 100. 

Correlation coefficient = -0.960. 



Figure S3. - Longitudinal velocity versus percent porosity for tungsten (ref. 2) with 
General Electric nominal 4 pm starting powder size. 

Velocity = -0.006 x percent porosity + 0.558. 

Percent theoretical velocity = -1 .13 x percent porosity + 100. 

Correlation coefficient = -0.91 0. 



Figure 54. - Longitudinal velocity versus percent porosity for tungsten (ref. 2) with 
General Electric nominal 18 pm starting powder size. 

Velocity = -0.008 x percent porosity + 0.554. 

Percent theoretical velocity = -1.38 x percent porosity + 100. 

Correlation coefficient = -0.992. 
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Figure 55. - Longitudinal velocity versus percent porosity for tungsten (ref. 2). 
Velocity - -0.006 x percent porosity + 0.533. 

Percent theoretical velocity = -1 .1 1 x percent porosity + 100. 

Correlation coefficient * -0.916. 



Figure 56. - Longitudinal velocity versus percent porosity for tungsten (ref. 2). 



Figure 57. - Longitudinal velocity versus percent porosity for UO2 (ref. 40). 
Velocity » -0.008 x percent porosity + 0.550. 

Percent theoretical velocity - -1 .49 x percent porosity + 100. 

Correlation coefficient - -0.997. 


44 


PERCENT THEORETICAL PERCENT THEORETICAL 

VELOCITY VELOCITY 






PERCENT POROSrTY 


Figure 60. - Shear velocity versus percent porosity for YBa 2 Cu 3 0 7 _ x (refs. 1 1 , 42 
to 45). 

Velocity - -0.002 x percent porosity + 0.313. 

Percent theoretical velocity = -0.766 x percent porosity + 100. 

Correlation coefficient = -0.814. 
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Figure 61 . - Ultrasonic velocity versus percent porosity for YBagCug 07 H( (refs 1 1 , 
42 to 45). 
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